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(57) ABSTRACT 

A method for sampling a signal for signal processing, such 
as calculating the coefficients of a Fourier or other transform 
of the signal. The sampling occurs at sampling points which 
are the union of sets of points, each set being points 
separated by regular intervals of l/p n where p„ is a prime 
number. Where the signal is not accurately reconstructed 
from the sampled values the method recurses and increases 
the set of distinct primes. The invention relies upon the 
incommensurate nature of different primes to prevent the 
invention from degenerating into an inefficient sampling 
technique such as using regular intervals would provide. The 
scale, termed a D scale, is applicable to problems in diverse 
domains such as signal processing, including digital signal 
processing (DSP), image processing, scientific and engineer- 
ing computational applications, data acquisition and statis- 
tical data analysis. Algorithms are provided that process an 
N dimensional Fourier and other kernel transforms of a 
multi-dimensional function using the non-uniform samples 
of the D scale along with Fast Fourier Transforms. 
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DIAGRAM 6: 
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DIAGRAM 1: 

UNIT INTERVAL, 1D, SINGLE-RESOLUTION PRIME NUMBER 3, DRUCK SCALE™: 
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DIAGRAM 2: 

UNBOUNDED, 1D, SINGLE RESOLUTION PRIME NUMBER 3, DRUCK SCALE 
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DIAGRAM 3: 
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DIAGRAM 4: 
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DIAGRAM 6: 

UNBOUNDED ONE (1) DIMENSIONAL MULTI-RESOLUTION DRUCK SCALE™: 
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DIAGRAM 7: 

ONE (1) DIMENSIONAL MULTI-RESOLUTION DRUCK SCALE™ BOUNDED BY 
SAMPLING TIME: 
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« CALCULUS FOURIERTRANSFORM // MATHEMATICA SPECIFIC. 
PULSEftJ = EXP[-.5t] UNITSTEP[t] // NOTE: e--5t FOR t >=0, 0 FOR t<0 

PLOT[EXP[-.5t] UNITSTEP[t], {t, -1 , 15}] // NOTE: PLOT FROM t=-1 TO t=15 
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SPECTRUM = FOURIERTRANSFORM [ PULSE[ t ], t, w ] 

1 // NOTE: w=2 n f WHERE f IS FREQUENCY. 

0.5 - 1 w // NOTE: CALCULATED SPECTRUM WHERE I IS THE IMAGINARY 'i' 

PLOT[ABS[SPECTRUM], {w, -10, 10}] 

// NOTE: CALCULATED SPECTRUM'S AMPLITUDE FROM |w|<10 
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- GRAPHICS - // NOTE: EXACT AMPLITUDE SPECTRUM 



FIG. 10 
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MATHEMATICAL BUILT-IN "FOURIER" FUNCTION PERFORMS JUST 
THIS CALCULATION: 

DATA40 = TABLE[PULSE[0.025n], {n, -10, 50}]; 

// NOTE: UNIFORM SAMPLING OF SIGNAL AT THE RATE OF 40/SEC 
(OR 1/40 INTERVAL). 

LISTPLOT[ABS[FOURIER[DATA40] ], PLOTJOINED-> TRUE, 
PLOTRANGE->{{0.,10.},{0. 1 2.5}}] 

// NOTE: AMPLITUDE SPECTRUM. 
" 2.5 r 

2 i- 

1.5 : 
1 \ 
■ 0.5 

2 4 6 8 10 
- GRAPHICS - 

FIG. 11a 




II NOTE: SAME AMPLITUDE SPECTRUM AS A FUNCTION 
OF co, ACROSS ENTIRE RANGE. 
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RESULTANT FREQUENCY D SCALE SUBSCALES - AMPLITUDES 
CLEAR[ g ] 

DO[g[m] = N[f[m]] 1 {m, LENGTH [p]}]; 

« GRAPHICS 'MULTIPLELISTPLOT' // NOT: MATHEMATICA SPECIFIC. 

LISTPLOT[ g [ 1 J, PLOT JOINED -> TRUE, PLOTRANGE -> ALL] 

2r 



1.5 



1 - 




- GRAPHICS - // NOTE: AMPLITUDE PLOT OF SPECTRUM 
CONTRIBUTION FROM D SCALE PRIME 2. 

FIG. 12 

LISTPLOT[ g [2], PLOTJOINED-> TRUE, PLOTRANGE-> ALL] 

2r 




5 10 15 20 25 30 35 

- GRAPHICS - //.NOTE: AMPLITUDE PLOT OF SPECTRUM 
CONTRIBUTION FROM D SCALE PRIME 5. 

FIG. 13 
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LISTPLOT[ g [ 3 ], PLOTJOINED-> TRUE, PLOTRANGE-> ALL] 
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MULTIPLELISTPLOT[ g [ 1 ], g [ 2 ], g [ 3 ], PLOTJQINED -> TRUE, 
PLOTRANGE->{{0.,10.},{0.,2.5}}] 
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MEASUREMENT SCALE FOR NON- 
UNIFORM DATA SAMPLING IN N 
DIMENSIONS 

This application claims the benefit of provisional appli- 

* cation No. 60/U5/752 filed Jan. 13, 1999. 

FIELD OF THE INVENTION 

This invention relates to methods for sampling objects. In 
particular it provides a system implementing a rigorous 
mathematical framework and provably correct practical 
algorithmic implementations for digitally recording a signal 
having any number of degrees of freedom. 

BACKGROUND OF THE INVENTION 

A physical or virtual object is recorded for application 

* specific analysis by sampling or measuring it at various 
points. An object might be a one dimensional electrical 
signal, price history or other time series, a two dimensional 
image, a three dimensional vector field, or an N dimensional 
surface in, say, configuration space, momentum space, or 
phase space, all useful in physics. One important example is 
the current use of uniformly spaced samples of signals or 
images to analyze these objects for subsequent processing 
like filtering or convolution. This analysis involves calcu- 
lating approximations to the Fourier transform, wavelet 
transform or other related transforms. The aggregate of 
sample points generated by the measuring process creates a 
mosaic from the original physical object. The degree of 
congruence of the mosaic to the original depends on the 
number of sample points used as well as on the judicious 
location of the sample points. 

A classical example of such is recording of an audio signal 
for replay by sampling the signal amplitude at uniform time 
intervals less than half the period of the highest frequency 
present in the signal. The so-called Nyquist theorem assures 
that such a recording enables reconstruction of the original 
signal. If the uniform sampling occurred at a lesser rate, in 
many cases an attempt to replay the sampled signal would 
. result in a distortion termed "artif acting* * in which harmon- 
ics of the signal would not correctly reproduce. 

Limited resources place constraints upon the generation 
of sample points. Therefore sampling is often a fairly 
expensive operation in terms of system resources consumed. 
Depending on the application, the number of sample points 
generated might be constrained by the resolution of the 
measuring equipment used (cost constraint), by the amount 
of storage available to save the samples (space constraint) as 
well as by the observation time window of opportunity to 
obtain the samples (time constraint). Thus, a set of samples 
can be regardedas a somewhat scarce resource. Therefore, 

* each sample point should have impact and provide a maxi- 
mal contribution of information content to the whole 
mosaic. 

Since there is a cost associated with the generation of each 
sample point, some optimization of sampling is generally 
required. Particularly in the case of objects having periodi- 
cally recurring patterns there is much redundancy in samples 
taken in uniform intervals. Or, more often, an object might 
have mostly unchanging features with little useful informa- 
tion content. But, say 10% of it contains volatile features of 
great interest. This has led to the use of non-uniformly 
, spaced, or irregularly patterned, data samples. 

For instance, a signal might exhibit rapid change (or 
volatility) in one concentrated region, but otherwise be fairly 
smooth. Or, an image might exhibit sharp edges as well as 
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regions of unchanging, or predictable, patterns. In both 
cases, many sample points are required in certain concen- 
trated intervals or areas where much information can be 
gleaned from the complex changes occurring there, while 

S few sample points are needed elsewhere where little useful 
information can be extracted. 

The number of sample points required using the standard 
uniformly spaced samples rises dramatically as the number 
of dimensions increases. Thus if a one dimensional signal 

30 requires O(N) samples, then a two dimensional image would 
need 0(N 2 ) samples, a three dimensional 0(N 3 ) 3 samples 
and so on. (A quantity Q(x) is said to be of order N*, written 
0(N*), if Q(x)/x* is bounded. 
What is required is a method to specify the manner of 

j5 sampling an object wherein each the totality of samples 
provides comprehensive and non-redundant information 
about the object, while minimizing the number of samples 
required, especially as the dimension of the sampled object 
increases. 

2Q The following methods are generally known in the prior 
art: 

Accept the limitations of uniform sampling and utilize 
Fourier Transform approximation techniques like the 
Discrete Fourier transform and its FFT implementa- 

25 tion ; . 

Approximate, generally via interpolation, non-uniformly 
spaced samples by uniformly spaced samples of suffi- 
ciently high resolution in critical regions tailored to the 
particular object; 
30 Using pseudo-random numbers, generate an irregular 
scale for measurements. A Monte Carlo method is 
typically used to implement this; 
Use multi-resolution wavelets; and 
Adaptive ly sample data. 
35 The Fourier Transform plays a dual role in N dimensional 
data recording. First consider a one dimensional signal 
representing time variations of some variable of interest, say 
voltage or price. The Fourier expansion of a signal is capable 
of representing an analog signal by a sequence of complex 
40 numbered Fourier coefficients, the amplitudes and phases of 
the Fourier components, whereby the content of a signal 
with a bounded frequency spectrum is accurately repre- 
sented by a discrete data series. This itself results in effective 
data compression, especially for signal having significant 
45 periodic components. Secondly, the determination of the 
Fourier coefficients of the signal is determined by sampling 
of the signal, which brings one back to the need for the 
present invention. 

The standard technique is to employ a Discrete Fourier 
50 Transform (DFI) as a substitute for the Fourier Transform. 
Furthermore the FFT or Fast Fourier Transform has been 
developed, which is an efficient implementation of the DFT. 
The DFT, therefore is an approximation that inherently 
distorts the spectrum generated by the Fourier transform. 
55 Nevertheless, it is tolerated as the most commonly used 
digital and image processing technique because it is ana- 
lytically tractable and readily implemented on diverse 
hardware/software platforms. 
The approximation, generally via interpolation, of non- 
60 uniformly spaced samples by uniformly spaced samples of 
sufficiently high resolution also distorts the calculated spec- 
trum. But it is also potentially wasteful of sample points 
since it requires uniformly spaced samples as anchors 
around which irregularly spaced samples are interpolated. It 
65 is also difficult to implement its mathematics. Indeed, the 
overwhelming general practice is still to use the first 
approach. 
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The use of pseudo-random numbers is analytically or 
mathematically difficult to describe, if not pragmatically 
intractable. That is, it is not readily describable 
mathematically, say as a set of equations. A mathematical 
description is generally needed in practice, to develop algo- 5 
rithms that automatically map sample points to the irregu- 
larly spaced random points. It is also computationally com- 
plex because the, deltas, or differences between each sample 
point and neighboring random points, generally must be 
recalculated because the random points used might not be 10 
reproducible. 

Another, relatively new, general technique is multi- 
resolution analysis. This involves splitting the stream of 
possibly non-uniformly spaced data into streams of different 
resolution. In some circumstances, this technique is also 15 
referred to as "decimation". Wavelets also use a multi- 
resolution analysis. The use of multi-resolution wavelets 
tries to solve the problem of mathematically decomposing 
non-stationary or highly irregular signals and images, into 
simpler components on multi- tiered scales. But the approach 20 
is generally difficult to implement in practice. Although 
there have been claims for very specific applications, its 
effectiveness vis-a-via standard DFT based techniques for 
general applications is not apparent. Indeed, the 1965 
Cooley-TUkey FFT algorithm immediately paved the way 25 
for the revolution in medical imaging that accelerated in the 
1 970s. Such a similar impact is not apparent, even more than 
13 years after the framework for wavelets was articulated in 
a coherent manner. 

There were other attempts preceding wavelets, at calcu- 30 
lating the Fourier transform of non-stationary objects, nota- 
bly the Gabor Transform. But again, these techniques were 
not general enough to be of use off-the-shelf, as the FFT. 
Indeed, a major motivation underlying wavelets was to 
supplant the Gabor transform method with a multi- 35 
resolution approach. 

Adaptive sampling has been used in many specialized 
applications in connection with the FFT. Its use of special- 
ized and sometimes customized sampling algorithms (to 
decide on how to adapt to data changes) limits its usefulness 40 
as a general purpose tool. It requires much a priori knowl- 
edge of the specific data set being sampled. 

U.S. Pat. No.. 5,708,432, appears to disclose a sampling 
digitizer with a sampling master clock divided by a factor k 
which is derived from prime integers. The '432 patent is 45 
concerned with sampling a periodic function, where the 
number of samples desired during each cycle exceeds the 
response ability of the sampling system. A technique is used 
in which fewer samples are taken from successive cycles of 
the signal and are combined in order to simulate the effect 50 
of many samples taken during one cycle. To do this, a clock 
is divided by an integer K which is varied according to an 
algorithm. N consecutive samples are taken over M cycles 
of the signal. Coherent signaling is achieved when the signal 
frequency and the divided clock frequency are relatively 55 
prime. The particular algorithm finds the relatively prime 
fraction by selection from a Farey series of a given order. 
Unlike the present invention, the '432 invention selects a set 
of samples separated by constant interval. The use of rela- 
tively prime ratios to guarantee that the same points are not 60 
selected on a perfectly periodic signal is different from the 
present technique, which would work to sample even a 
pulse, and does not require that the signal be periodic. 

An abstract in the prior art makes reference to sampling 
employing prime numbers. Proceedings the European Signal 65 
Processing Conference held in Trieste on Sep. 10-13, 1996. 
The abstract states: 



We demonstrate that it can be advantageous from a 
computational point of view to use a two-stage real- 
ization instead of a single-stage realization for sample 
rate conversions with prime numbers. One of the stages 
performs a conversion by a factor of two .... The other 
stage changes the sample rate by the rational factor N/2 

U.S. Pat. No. 5,243,343 uses variable time interval data 
sampling. See col. 5, line 18-19. U.S. Pat. No. 5,444,459 is 
essentially identical with the '343 patent. 

Other prior art shows 

1. non-linear sampling without reference to the prime inte- 
ger algorithm of the present invention (U.S. Pat. Nos. 
4,188,583; 4,142,146; 4,903,021;and 5,229,668 [FIG. 3]) 

2. the use of prime numbers in Hartley transform processors 
(U.S. Pat. No. 4,062,060; Boussakta et al., 24 Elect. Ltrs. 
No. 15, pp. 926-28 (Jul. 21, 1988); 136 I EE Proc. Pt. G, 
No. 5 (October 1989) pp. 269-77.) 

3. Other patents of interest include U.S. Pat. No. 5,712,809 
entitled "Method and apparatus for performing fast 
reduced coefficient discrete cosine transform"; U.S. Pat. 
No. 5,682,524 entitled "Databank system with method for 
efficiently storing non-uniform data records"; U.S. Pat, 
No. 5,497,152 entitled "Digital to digital conversion using 
non-uniform sample rates"; U.S. Pat. No. 4,271,500 
entitled "Device for converting anon -uniformly sampled 
signal with short-time spectrum to a uniformly sampled 
signal"; U.S. Pat. No. 5,712,635 "Digital-to-analog con- 
version using non-uniform sample rates*';. U.S. Pat. No. 
4,999,799 entitled "Signal processing apparatus for gen- 
erating a Fourier transform"; and U.S. Pat. No. 4,969,700 
entitled "Computer aided holography". 

4. 

BRIEF DESCRIPTION OF THE INVENTION 

The invention relates to signal processing, primarily digi- 
tal signal processing employing N dimensional data sam- 
pling required to extract information about the object for 
subsequent processing. In order to extract information from 
a signal it is sampled at discrete points and the sampled 
values are fed into calculation algorithms. For example the 
sampled values may be used to calculate the coefficients of 
a Fourier transform of the signal by known algorithmic 
techniques, thereby extracting useful information about the 
Fourier spectrum of the signal. Other transforms may be 
calculated using other kernel functions to generate Radon, 
Hartley, or discrete sine and cosine transforms. The discrete 
points at which sampling occurs may for example be dis- 
crete times for a signal that varies with time. As a specific 
example of the utility of the invention, it is utilized to 
process an N dimensional Fourier and related transform 
using the non-uniform samples of the D scale and FFTs. This 
decomposition of a Fourier transform into a series of FFTs 
will be detailed later. 

This invention is concerned with the selection of the 
values of the discrete points at which sampling occurs. Since 
most signals are composed of a few periodic components, 
there is much redundancy when the sampling takes place at 
regular intervals. Sampling at regular intervals also leads to 
artif acting, i.e. false data caused by the confusion of an 
interval with its harmonics. On the other hand it is algorith- 
mically simpler to sample at regular intervals and it is 
necessary to keep the calculations manageable. 

This invention provides a method including algorithms to 
generate an irregularly spaced set of sampling points which 
still provides the practical calculation advantage of uniform 
sampling. This methodology can also be tuned or configured 
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to generate an approximation of a truly (pseudo)random set familiar ruler scale as a point of reference. The new scale is 
of points. Further, it can be tailored such that the resolution based, fundamentally, on the existing mathematics of set 
or width between neighboring points in the set is always theory, real analysis of the continuum, and number theory, 
within pre -specified maximum and minimum bounds. That That is, it applies theoretical knowledge from set theory, real 
is useful when sampling different data sets, each with a 5 analysis and number theory to solve practical problems 
different number of significant digits. For a one dimensional related to data sampling. The theoretical framework pro- 
situation, i.e. where the signal is a function of one parameter vides the necessary rigorous foundation to guarantee that 
such as time, the set of sampling points is composed as the this new scale yields provably correct results in all applica- 
union of n sets of points, each set being points separated by tions. 

regular intervals of l/p„ where p n is a prime number. Each The new scale, herein termed a D scale, is applicable to 

set is denoted a "scale" for prime number p„. The union is problems in signal processing, including digital signal pro- 

therefore characterized by a set {P n } of distinct prime cessing (DSP), image processing, scientific and engineering 

numbers. The union is intrinsically non-periodic. And for computational applications, data acquisition and statistical 

even a small set of distinct primes the distribution of points data analysis, to mention a few application domains, 

in the union appears almost random. This union is used to The D scale can also be configured or tailored to match 

determine where the signal is sampled. The points at which the varying bandwidth, resolution or accuracy, and point 

sampling takes place is not dependent on the. values of the distribution requirements of the widely disparate application 

signal If the signal cannot be accurately reconstructed from domains mentioned above. Specifically, the minimum and 

the sampled values one merely recurses and increases the set maximum resolutions of any particular D scale can be 

of distinct primes, which is always possible because of their on guaranteed by an algorithm discussed later, that determines 

infinite number. Alternately, there might be a-priori knowl- 20 ^ P nme ^ ales 1 to mch ! dc when a * cale * ™™t™ted. 

edge of an object's information bandwidth. Then the prime Further ' * c f *f t °™ * constructed with prime scale 

& , . J 4 ... * r • * i_ components such that the distribution of its points is made 

numbers used to generate toe set of points can be pre- somewhat random for certain applications . 0r , the points can 

determined, such that ttie highest prime number just exceeds be distributed t0 accentuate the multi-frequency components 

the known, effective, bandwidth. For example audio apph- 25 of m application domain> suiting in more points towards 

cations generally require bandwidth ranges of 5-20 KHz lne center of a unit ^^1, and tapering off at the edges. 

Therefore there is no need to use pnme numbers of order 10 ^h^ flexibility in generating point distributions is effected 

when this application domain only uses an effective band- by including more or fewer prime scales in a D scale, 

width of order 10 . j^c invention employs a method, or an algorithm, to 

The invention may be implemented in a real-time mea- 30 perform nonuniform sampling. It can be embodied or imple- 

surement operation by running a series of clocks ticking at mented as a product (hardware based, software based or 

intervals characterized by l/p„ for the n r/l clock, and sam- combined hardware/software) to implement nonuniform 

pling the signal at each clock tick. For a more than one sampling. 

dimensional situation, where the signal is characterized by possible mmmixM applications are the follow- 

multiple parameters, each parameter is treated as in the one 35 m g. 

dimensional case, and the sampling space is treated as the R ea l-time Data Acquisition 

Cartesian product of one in dimensional spaces. The inven- j^e D scale ^ along with its samp ii ng algorithm, improves 
tion relies upon the incommensurate nature of different the efficiency of real-time data acquisition. Data acquisition 
primes to prevent the invention from degenerating into an fe at ^ heart of many processes m such diverse fields as 
inefficient sampling technique such as using regular inter- 40 seismology, oceanography, medical imaging and process 
vals would provide. control, to name a few. The D scale and its algorithm 
Although the invention is characterized by the use of samples judiciously, say when volatile data is detected, and 
primes, there are other techniques that may be used to not wnen the data is uninteresting. In contrast with standard 
generate the sampling points. The key features are that the techniques that leads to reduced cost of data acquisition 
points are derived from the union of sets of points generated 45 smce f cwer sam ple points are generally used. The effective- 
by simple calculations and the union avoids a periodic ness 0 f lne sampling process is improved since useful 
structure. Another variation of this sampling technique is to samples are generated, not the many extraneous samples of 
offset the points associated with each prime scale. The a uniformly spaced procedure, 
onsets might be random or also prime numbers. The intent Post- acquisition Data Analysis 

would be to better capture phase information in data analy- 50 The D scale can also be used to analyze data already 

sis- generated by a data acquisition or sampling device. For 

An efficient method is presented for sampling and thereby example, the D scale is used to generate a closed form 

enabling the recording of, or measurements of objects, expression for the Fourier Transform (described later), 

which is successful regardless of the irregularity or period- Using the closed form expression, it is then possible to 

icity of the object. The sampling occurs at values of an 55 quantify the bias in not using points in regions that are not 

independent variable. For example, where the sampled varying enough to warrant attention. Therefore, it can effec- 

object is an electrical signal, the independent variable could tively filter the acquired data or samples to focus on the more 

be time. Wherethe sampled object is a multidimensional significant data sections, thereby reducing the quantity of 

physical object, the independent variable could be a spatial data that must be processed to characterize the object 

coordinate parametrizing the space in which the physical 60 sampled. In this manner the D scale acts as a means for data 

object is embedded. Collectively, the points or values of the compression. This has application in the analysis of large 

independent variable at which samples are taken is termed a quantities of data produced by NMR or tomography devices, 

"scale". This terminology is consistent with the usual asso- and reduces the burden of analyzing such data. The same is 

ciation of a scale with the ruler used to measure quantities. true of the large quantity of data produced by seismography 

This new measurement scale is presented with special 65 measuring apparatus, as an example. Further, other applica- 

properties that makes it particularly effective for sampling tion domains produce resonance shapes which might contain 

non-uniform data. The new scale is contrasted with the only 10% useful information content. 
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Transforming Irregular Data into Periodic Components 

The D scale can be used to partition an original, irregu- 
larly patterned data set to create multiple components. Each 
component contains points which, although not periodically 
patterned, are all reciprocal multiples of the same prime s 
number 1/p^, which is what bounds them to the same prime 
scale. Then, because the points on each prime scale is 
bandlimited by p*, the set of points on each prime scale can 
be decomposed via standard Fourier analysis into periodic 
components. Thus, the D scale enables a two stage decom- aQ 
position. The first stage decomposes the original set of 
irregularly spaced points into sets of uniformly spaced and 
bandlimited points, each associated with a prime scale. Then 
the second stage decomposes each bandlimited decomposi- 
tion into the usual periodic components of a Fourier series. 
Also, each partition has a different resolution (l/p^such that 15 
the components form a nested set of resolutions. 
Characterization of Noisy data 

The D scale can filter random noise into periodic com- 
ponents thereby providing a useful characterization of the 
data, which reduces the effect of random noise in the data. 20 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 shows a bounded, ID, single resolution prime 
number 3, D scale. 

FIG. 2 shows an unbounded, ID, single resolution prime 25 
number 3, D scale. 

FIG. 3 shows a unit interval, ID, multi-resolution D scale 
comprising a composite scale using the prime numbers 2, 3 
and 5. 

FIG. 4 shows an unbounded, ID, multi-resolution com- 3Q 
posite D scale. 

FIG. 5 depicts the steps of the process for calculating the 
Fourier transform spectrum of an object. 

FIG. 6 depicts the 1 dimensional component of an N 
dimensional D Scale comprising potentially an unlimited 35 
number of prime resolution scales. The union of points on 
these scales creates an irregular pattern of points on the 
composite D Scale. 

FIG. 7 depicts a one dimensional multi-resolution D scale 
bounded by sampling time. The spectrum's boundary is 
determined by the prime number of the resolution used. FIG. 40 
7 shows this effect. 

FIG. 8 depicts the reciprocity theorem. 

FIG. 9 depicts a benchmark signal used to evaluate the 
efficiency of the invention. 

FIG. 10 depicts the exact spectrum of the benchmark 
signal of FIG. 9. 

FIG. 11a depicts the amplitude spectrum of the bench- 
mark signal arrived at by using Mathematical Fourier 
function. " 50 

FIG. lib depicts the amplitude spectrum as a function of 
co, across its entire range. 

FIG. 12 depicts the amplitude plot of the spectrum con- 
tribution from the D scale using prime 2. 

FIG. 13 depicts the amplitude plot of the spectrum con- 55 
tribution from the D scale using prime 5. 

FIG. 14 depicts the amplitude plot of the spectrum con- 
tribution from the D scale using prime 11. 

FIG. 15 depicts the spectrum from D scale primes 11, 5, 
and 2 contributions. 60 

FIG. 16 depicts the spectrum from D scale primes 29, 19, 
and 11 contributions. 

DETAILED DESCRIPTION OF PREFERRED 

EMBODIMENTS 65 

The D measurement scale has three components; a set of 
scale points used-to select data samples, a set of pre- 
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calculated distances between nearest neighbors and an algo- 
rithm that efficiently associates or maps sample points to the 
scale points. The rigorous mathematical underpinning of the 
D scale involves two fundamental mathematical concepts; 
the notion of a one-to-one correspondence between points in 
two sets and that of the measure of a set The measure of a 
set corresponds roughly to the inherent concept of size, e.g. 
length or area, but is defined precisely to cover sets whose 
inherent size is not obvious. The Lebesgue measure of a one 
dimensional set is a familiar example of the use of the 
concept of measure, and permits assigning a measure for 
such oddities as the set of all irrational points, etc. 

The familiar ruler uses a subset of rational numbers (i.e., 
numbers that can be expressed as fractions, typically mul- 
tiples of some high inverse power of two, which defines the 
resolution of the scale) for measurements. The rulers sam- 
pling algorithm finds the multiple of the inverse power of 
two that is nearest to the sample point. For example, an 
irrational number is approximated by a rational number of 
the required resolution. 

The D scale satisfies several criteria for the utility of a 
measurement scale. It has 

a. A degree of regularity, or of a uniform pattern, of the 
scale points for ease -of implementation, that is to say 
the definition of the scale points is simple to calculate, 

b. A degree of irregular pattern of scale points to accom- 
modate clustering of data points, that is to say the scale 
points exhibit a highly non-uniform distribution, 

c. Multi-resolution capability to accommodate different 
orders of magnitudes in sample point resolutions, 

d. Minimal number of scale points needed for effective 
measurement. 

Note that the familiar ruler meets criteria a and c, pro- 
viding uniformly spaced markings and multi-resolution as 
powers of 10. In contrast, a scale whose points are randomly 
selected only meets criteria b above. The D Scale meets all 
four criteria indicated above, therefore making it most 
effective for nonuniform data sampling. 
Working Description of D Scale Invention 

The D measurement scale is described both analytically 
and geometrically. Then, rigorous mathematical proofs are 
presented on the correctness of its operation. After that, the 
scale's sampling technique is described and its utility 
proven. Finally, useful and commercially viable examples 
are provided illustrating uses of the D scale. 

The N-Dimensional D scale consists of selected points on 
the N dimensional rational number coordinate space (or 
grid). This scale will be shown to meet each of the criteria 
set forth for the D scale. It will also be shown to be effective 
in sampling and/or processing non-uniformly spaced data. 

The set of measurement points in a D scale is described 
in stages. First, the one dimension (ID) scale is described. 
This is then generalized to a description of the N dimen- 
sional scale. 

The ID measurement scale is multi-resolutional, i.e. it is 
comprised of multiple scales of different resolutions. The 
ID, unit interval scale with a single resolution is first 
described. Then, the unbounded 1 D scale with a single 
resolution is described. That is followed by the multi- 
resolution, unbounded, ID scale. 

Consider a one dimensional bounded single resolution 
scale, e.g. the set of points on the ID, half -closed interval [0, 
1), that are multiples of some prime number. Or more 
formally, consider the set of points S^ ^-{x € [0, 1): x-o/p, 
n=l..p-l, pe {prime numbers}}, which is depicted in FIG. 
1 for the prime number 3. 
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Next, extend this definition to an unbounded set of points I have accomplished rigorous proofs of three theorems 
by replicating the set of points on [0,1) across the entire that demonstrate properties establishing the utility of the D 
number line. More formally, consider the set of points S p ={x scale. The first property is the scale points' uniqueness. The 
e (-oo 5 oo):x=j+n/p, j=-oo . . , oo, n~l..p-l, p e {prime second is the scale points' ability to completely cover a line 
numbers}}. S p is referred to as the p rA prime resolution D S (in ldimension) or region (in N dimensions) without over- 
Scale. This is depicted in FIG. 2. lap. And the third property is the scale points' capability to 

Next consider a multi-resolution bounded D scale. Con- provide any desired resolution through nesting of prime 

sider N ID, unit interval, prime resolution scales, each resolution scales. These properties are direct consequences 

corresponding to a different prime number. A multi- of the intrinsic properties of prime numbers. Indeed, that is 
resolution D Scale is then the superposition or union of 10 why only prime numbers are included in the D Scale. But 

measurement points from each of these N prime resolution more generally, any other set of numbers should have these 

scales. Formally, using the previous definition of S^ 0 the characteristics as well. Specifically, using random or prime 

ID multi-resolution D Scale is: S^ 0 ^-{xe US^ 0 ^ip-p, . . . numbered offsets should be sure to obey these theorems. 

Pn> P/ € {prime numbers)), or equivalently, S^^-US^^ These theorems underlie the validity of the scale and provide 

where p-p, . . . p M , p ( - c (prime numbers)). This is shown in 15 the necessary confidence that the scale always performs its 

FIG. 3. measurement function correctly. 

Next, consider N ID, unbounded, prime resolution scales, Theorem 1 is the Fundamental Uniqueness Theorem 

each corresponding to a different prime number. A multi- Each D Scale point is in one and only one prime resolu- 

resolution D Scale is then the superposition or union of tion scale. That is, for any pair of N-th dimensional prime 

measurement points from each of these N prime resolution 20 resolution scales S ( ^ and S, in the D Scale it follows that 

scales. Formally, if S p is as defined previously, then the ID S it u O S yA =0 where I, j e {prime numbers}, 

multi-resolution D Scale is S,«{x e US^p-p, . . . p„, p, e Theorem 2 (One Dimension) 

(prime numbers }), or equivalently, S^US^, where p=p, ... The set of D Scale point neighborhoods covers the 

P*> Pi e (prime numbers }). This is shown geometrically in number line (-o°,<») (i.e., Rl), And the covering does not 

FIG. 4. 25 contain any overlapping neighborhoods: UNp (V =(-o°,o°). 

Summarizing, each measurement point in the multi- Corollary 1 (1 Dimension): Any sample point can be 

resolution ID D Scale is located by its three coordinates assigned or mapped to one and only one D Scale point. 

(p,ij), where p is the prime resolution scale in which the Theorem 2 (N Dimensional) 

point is located, i is the integral interval on the prime The set of D Scale N dimensional point neighborhoods 
resolution scale p, and j is the fraction of the prime number 30 covers the R^ continuum. And the covering does not contain 
in that i** interval. Therefore, if a measurement point's any overlapping neighborhoods. U^^co^. (These neigh- 
location on the rational number based ruler is d, then the borhoods are defined below.) 

relation of these scales is d=i+j/p. Corollary 1 (N Dimensions): Each sample point can be 

Progressing now to multidimensional examples, a single assigned to one and only one point neighborhood in the D to 

resolution, n dimensional, D scale is defined for the bounded 35 Scale. 

unithypercube. First, consider Cartesian cross product of the Definition: The resolution (or span) of a D Scale is the 

measurement points on the one dimensional half-closed maximal distance between adjacent points such that they are 

interval [0, 1), namely, S pJv j 0j S^^oS^o ^x . . . x distinguishable. The gaps between adjacent D Scale points 

S/jfo.i) (N times). An N dimensional point on the unit varies. But there are applications where the maximal gap 

hypercube D Scale is then described as X e IT S p ^ 0 3) -»X= 40 between adjacent points must be guaranteed to within a 

(x ; . . . Xjy)-*Xy € S^q fl y specified resolution. 

This has also been extended to the unbounded hypercube. Theorem 3 (Infinitely Scalable Resolution) 

First, consider the measurement points on the one dimen- The span of the D Scale decreases to 0 as the prime 

sional half -closed interval [0, 1), previously defined as S^ 0 resolution tends to oo. Conversely, the prime resolution tends 

i). Define the Cartesian cross-product of N of these 1 45 to oo as the span decreases to 0. This theorem guarantees that 

dimensional unbounded scales as follows: S^U S p «S p the resolution of a D Scale can be adjusted to as small a 

x . . . xS^ (N times). An N dimensional point on the value as required, by incorporating higher prime resolution 

unbounded unit hypercube D Scale is then described as X e scales. 

S^-*X ellS^ X=(x; . . . Xj^-^Xy e S^. Corollary 4: A prime resolution scale can always be found 

This in turn is expandable into a multi-resolution, N so to implement a specified maximum gap on an N dimensional 

dimensional, D Scale. For example, an N-dimensional D scale. Additionally, a theorem is proved Theorem 4: A 

multi-resolution scale on a unit hypercube is created by the necessary condition for neighboring scale points to always 

union of N dimensional, unit hypercube, prime resolution have a minimum interval, e, between them, is that: 
scales, for each prime number included in the scale. 

Formally, the N dimensional, multi -resolution, D Scale is 55 \Vpr^lPui\<^ 

S M o,i)-{X e US^ 0(1) :p- P/ . . . p„, p,. e {prime numbers}}, f 

or equivalently, S M0 n=US^ 0 n where p=p, . . . p„, p, e tor J? 1 . 1 p < P nm e fcaJes. - , . ^ . 

{prime numbcS}}Tiquivffiil?y > S m J<l { . . . ^K'u ^ * ***** pr ° Ven USing mathematical inducUon > 

S^q 3) where p=*p, . . . p„, p, e {prime numbers}}. since. 

An N-dimensional multi-resolution scale on an 60 |i//^VaI<<H"/a>-™/PiI<«<Po, «<*v 
unbounded hypercube is created by the union of N 

dimensional, unbounded, prime resolution scales, for each If 
prime number included in the scale. Formally, the N 

dimensional, multi-resolution, D scale is: S^{X e US^p- \VPrVPM\<t-~\*fp k -'*lp\<ui<p k , m< pl 

P/ • • • P«> P f - € {prime numbers}}, or equivalently 65 
S^USpj^p^p, . . . p„, p t e {prime numbers}, or equivalently, 

Sjv-{(x, • • . XAr):x y c US^p-p, . . . p„, p,- e {prime numbers}. \VPi+i-Vp<J<^\VPi+i-Vp l + 2 \<*-*\n/Pk-m/Pt\ 

<€/i<p t} m<p t 
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by the induction condition. And 

i+V m< Pi 

Therefore In/p^-m/p^e n<p jfc4 . 1 , m<p, l<k+l. Therefore 
the induction hypothesis is true for case k+1. 

Therefore a D scale can be constructed following an 
induction argument. Start with the lowest primes that satisfy 
the constraint. Then add a third prime that satisfies the 
constraint with the second prime, and onward. This con- 
structs a D scale layer by layer, or prime subscale by 
subscale. 

DETAILED DESCRIPTION OF A PREFERRED 
EMBODIMENT OF THE ALGORITHM FOR 
IMPLEMENTING A D SCALE 

Non-Uniform Sampling of Data Using the D Scale sam- 
pling can be performed in real-time, as each data point is 
individually received. Or sampling can be performed after a 
collection of data points is first gathered. In either case, an 
algorithm is presented to "measure" or sample data points 
using the D Scale. 

This algorithm is first presented for 1 dimensional data 
sets and then generalized to N dimensional samples. The 
algorithm is presented as a theorem, with the proof of its 
utility relegated to the examples. 
Theorem 5 (Sample Measurement Theorem in ID) 

A sample point x can always be assigned, or mapped, to 
one and only one point in a D Scale according to the 
following algorithm: 

Start from the lowest prime resolution scale in the D 
Scale. Determine if the sample point is in the neighborhood 
of any of the points in the D Scale. If it is, then assign that 
scale point to it. Otherwise, continue to the next prime 
resolution and perform the same search for nearest measure- 
ment point. The sample point must eventually find itself in 
the neighborhood of some point in the D Scale. Formally: 

for (P t - where P^p, . . . p„, p ( e {prime numbers in this D 
Scale}) 

Consider a sample point x, expressed as a decimal. 
Then x=a.b, where a is an integral part of the 
decimal, and b is the fractional part. The integral part 
of this sample, a, is the a-th interval on the prime 
resolution scale P,-. The fractional part, b, is used to 
determine which multiple of to use, as follows. 
Consider: (.b^P,- 

The floor function extracts the greatest integer j less 
than or equal to f: j=[(.b)*PiJ 

Therefore: j/P f =<b 

Also: .b=<G+l)/P t - 

Otherwise, j+l=[( b)*PiJ which contradicts the above. 
Thus: j/P,=<.b=<(j+l)/P 1 . 

The sample point's remaining distance from (j/P,) is: 

As verification, the sample point can be located on the a* 
interval of the ] th element of prime resolution scale P,-. 
Or: 

x=a+(j/Pl) +€ »=a. b 

Next, determine if the sample point is in the neighborhood 
of either i th or j+1 st measurement points of prime 
resolution ? t ; P Sj or P, v+1 respectively. 

Let S /V? =distance of P lV to its nearest right neighbor 
among the D Scale measurement points. 
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distance of P^ a to its nearest left neighbor among the 
D Scale measurement points, 
[f(e-<5,.y2) 

then measurement point P,-^ is selected as that sample's 
5 nearest approximation. That D Scale point is used in 

subsequent calculations to represent the original 
sample point. 
Else if (e-<V 1>t /2) 

then measurement point Py , is selected. 

10 

Else continue the same logic, using the next prime reso- 
lution scale. 

Theorem 4 (Sample Measurement Theorem in N 
Dimensions) 

An N dimensional sample point can always be assigned, 
or mapped, to one and only one N dimensional point in a D 
Scale according to the following algorithm: 

Start from the lowest prime resolution scale in the D 
Scale. Determine if the sample point is in the neighborhood 
of any measurement points in the D Scale by examining the 
ID neighborhoods of each of its coordinates, at the same 
resolution scale. If it is in the neighborhood of a scale point, 
then assign that scale point to it. Otherwise, continue to the 
next prime resolution and perform the same search for 
nearest measurement point. Hie sample point must eventu- 
25 ally find itself in the neighborhood of some point in the D 
Scale. 

Formally: 

For (P, where P^p,, . . . p„, p, e {prime numbers in this 
30 D Scale}) 

Consider an N dimensional sample point X, expressed 
as: 

X-{xl, X2, . . . Xjy) 

35 Then apply the ID algorithm to each coordinate x,- to 
that coordinate point's nearest ID scale point 
Coordinates at different dimensions x, will find their 

neighborhood at different prime resolutions x pk . 
Eventually each coordinate must find its neighborhood. 
40 The sample point will then be assigned to the scale point 
as: 

X-{xl, x2, . . . (x pl> x p2 . . . x pN ) 

General Properties of the D Scale 

45 FIGS. 6 and 7 summarize the key properties of the N 
dimensional D Scale, even as the 1 dimensional case is 
highlighted. One property is that the scale is composed of 
potentially an unlimited number of prime resolution scales 
each of which has points that are uniformly spaced. The 

50 union of points on these scales creates an irregular pattern of 
points on the composite D Scale . The availability of a 
countable infinite number of prime resolution scales leads to 
any desired resolution. Further, the scale can be tailored to 
the application domain of interest such as voice or images. 

55 Thus, the resolution between neighboring points can be 
made bounded, or the point pattern generated to be more or 
less random. 

Another property is that each scale point has a neighbor- 
hood about it that does not overlap with the neighborhood of 

60 any other scale point. This property leads to a measurement 
algorithm that is guaranteed to identify the neighborhood 
within which a sample point resides, in order to map it to a 
point on the scale. FIG. 6 illustrates this effect. Thus, it is not 
possible for a sample point to "fall-through the cracks" and 

65 not map to any scale point. 

Yet another property is that the spectrum of an object's 
(e.g., signal, image) sample points on any one prime scale is 
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necessarily band limited. The spectrum's boundary is deter- 
mined by the prime number of the resolution used. FIG. 7 
shows this effect. 

Contrasting this scale with the standard ruler, it requires 
far fewer measurement points than a ruler of the same s 
resolution or depth. Thus, at any depth measured by the 
number of decimal places (and therefore resolution), there 
are fewer ruler points per unit interval than with the standard 
ruler. 

Example of Implementation of the D Scale 3Q 

The D Scale can be completely implemented in software, 
in hardware or in both software and hardware. 

Consider the case of post-acquisition data analysis men- 
tioned earlier. Then, the scale, with its associated mapping, 
can be implemented, along with sampling, entirely in 
software, on already acquired data. Such is commonly the 35 
case in medical imaging or seismology, to touch on a few 
general application domains. Or, consider the case of real- 
time data acquisition, say in radar. Then it can be imple- 
mented entirely in hardware such that the sampling is 
automatically performed in real-time, as data is received by 20 
a data acquisition device. Or the scale's points can be 
mapped to the D scale in software, while its associated 
sampling algorithm is implemented in hardware. 
The example Includes: 

1. A set of prime numbers used to comprise the scale's 25 
multi-resolutions. The set is selected with the aim to bound 
the range of resolutions, according to an application domain 
constraint, with as few points as needed. 

2. Set of points in the scale as the union of points in each 
of the scale's resolutions. 30 

3. Distances between adjacent points on scale. It is saved 
in-memory or elsewhere. 

4. Sampling algorithm. 
Calculating the N Dimensional Fourier Transform Using the 
D Scale 

The underlying mathematical foundation of many diverse 
applications includes the Fourier Transform. These applica- 
tions are in such diverse domains as signal processing to 
include digital signal processing (DSP), image processing, 
scientific and engineering computational applications, data 
acquisition and statistical data analysis. But there are several 40 
difficulties in actually applying the Fourier Transform 
directly to measurement samples. 

Therefore in practice, several approximations of it are 
routinely used instead. They include the Discrete Fourier 
Transform which is typically calculated using the FFT 45 
family of algorithms, the Discrete Cosine Transform and 
others. There have been many attempts at calculating the 
Fourier transform of a function using irregular samples. 

The problem is then to calculate an accurate spectrum 
approximation of the following Fourier transform (using the 50 
time-frequency analog): 

given a non-uniformly spaced set of time samples. 

Part of the problem of calculating the transform is to 
provide feedback on how good the approximation is, given 
the set of sample points applied. The algorithm also can 
provide guidance on whether a redeployment of those 
sample points provide a more accurate spectrum. Further, it 
can be used to determine whether more sample points are 
needed to gain noticeably increased accuracy, or whether the 
point of diminishing returns has been reached. 

The present invention is used to provide a closed form set 
of equations, to calculate the Fourier transform of a complex 
valued function given a (finite or infinite) set of non- 
uniformly spaced sample points. The algorithm can also be 
applied without change to its procedural structure, to calcu- 



late similar transforms, such as, but not inclusively, the 
Harley transform, Discrete Cosine Transform (DCT), Fou- 
rier series. The algorithm is fast because it expresses the 
Fourier transform as a series of Fast Fourier Transforms 
(FFT). 

The Fourier transform is calculated by starting with first 
principles, from the definition of a Riemann Integral of a 
multi-dimensional function. The Riemann partition used to 
approximate the integral (say as an area under a curve), is 
determined by the set of measurement points in a time 
domain D scale. Also, the set of frequency components 
calculated is determined by a set of measurement points on 
a frequency domain D scale. Finally, the Shannon Sampling 
theorem is applied to calculate the approximate time domain 
function that corresponds to the just calculated frequency 
spectrum. 

The Fourier Transform of a function is expressed as a sum 
of band limited functions. Each of these band limited 
functions is derived from a subset of points in the D scale 
associated with one of its prime scales. Because of that, each 
of the bandlimited functions can be calculated efficiently 
using the standard FFT. Therefore, the Fourier Transform 
can be calculated with the same order of efficiency as FFT. 
Further, because each function component is bandlimited, 
each one can be calculated using the Shannon Sampling 
Theorem (which relies on the Nyquist Sampling Theorem to 
reconstruct an original object from its bandlimited sample 
points). Therefore, an approximation of the original object 
can be reconstructed from the set of all components, pro- 
viding the necessary feedback on the degree of likeness to 
the original object. Then, the sampling can be adjusted 
according to the desired accuracy, by changing the sample 
point deployment and possibly adding more samples. 

The steps of the process are depicted in FIG. 5 illustrating 
the iterative steps to approximate the Fourier Transform of 
an object (in N dimensions), using the sample points avail- 
able. 

This Fourier Transform calculation can only meaningfully 
calculate as many points in the frequency spectrum as are 
available sample points. That is, it is not possible to generate 
more information than what was provided. This constraint is 
readily expressed by the following matrix equation: 

Frequency vector=Fourier Decomposition Matrix 'Time Sample 
Vectorfwo . . . (i)jv_ t )=Afx[f 0 . . . t N ^\ 

This equation shows how the set of derived points in a 
spectrum is structurally related to the sampled point inputs. 
The detailed nature of this matrix structure follows. The 
points in the frequency vector are determined according to 
the Reciprocity theorem. This theorem states that: 

If A is the scale bound in one domain (typically time), Q 
is the scale bound in a second domain (typically frequency), 
and N is the total number of samples in the first domain, 
then: 

An alternate formulation is: 

If At is the sampling interval in one domain (typically 
time), Af is the interval in a second domain, (typically 
frequency), and N is the total number of samples in the first 
(time) domain, then 

Consider sampling a signal on prime scale p, for T p time 
intervals, where T p is a prime number. Then the boundary is 
65 T p . Then At-l/p where p is the subscale's prime sampling 
rate and the number of samples N«T p *p on subscale p. On 
the frequency side the roles are reversed, Af-l/Tp on sub- 
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scale T p . Let p represent the sampling rate of a D Scale's 
prime, subscale. 

This structure* satisfies the reciprocity theorem. 

At *A/-1/ J p" l/I>l/fT/p)-l/JV 

The set of derived frequencies is then related by the reci- 
procity theorem to the sample points by a one-to-one cor- 
respondence between time subscales and corresponding 
frequency subscales. 

These calculations then yield an expression for the Fou- 
rier Transform decomposition for each frequency subscale 
characterized by T^. for each corresponding frequency sub- 
scale: 

i r (2nm/T^>Real(2jm^>iTmagiiiary(2jtm/r^) 

or 

F(2mn/T pi )=\F(2jxmIT pl )| cxp (-t6) 

Where: 

lF(2jim/T^|=SORTlImagiiiaryC2jtm/7^) 2 +Real(23im/r />/ ) 2 ]- Ampli- 
tude 

e(2jtm/r p( )=Arctan(Imagmary(2^ 

IF: 

1 f*" 1 r T pnax 

F.iohn) = 2* J h{r)€~ "V dt = -n J o A(i)e""V dt 

THEN: 
\F{2xmlT Pi )\ = 

[m« prime) n ={ T pk VlH jT ( . 

t=(min prime) n=0 

where: 

The prime scales associated with prime numbers p ( - are 

components of a D Scale 
T p , is the observation length, in units of time, on prime 

scale p,.. 

T^ must be a prime number. T pmax is the maximum value. 

m=0, . . . T pi -1 for each frequency prime scale T pi 

where k=min to max prime scales. 

§ pkn %pk is the interval between a point on the D Scale 
located by (pj^n %p fc ) and its next, right, nearest neigh- 
bor on the same D Scale. (% is an operator extracting 
the non-integral part of a fraction.) 

{max prime} has a greater bandwidth than the function 
being transformed. Z' is a sum over the prime numbers 
in the scale. 

If the original signal, h(t)is expressed as a complex 
function h(t)=h^(t)+i h,<t) 
then: 

{max prime} n-(T pJt * Pjfc) -l 
Real(2rnn/T 16^ %/ A(o/p*) cos (2n{T/r>m/ 

Tp/p*)+ <V» ^/AOVp^siQ 
(2*(X/T>m/V PA ) 
k={min prime} n=0 
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Imaginary(2jmi/r^)=***2: 26^ %p M^Pid cos (2ji(T / 
T Jnm/T /p^+8^ ^/n/p^sin 

Discussion of Closed Form Solution 
S Consider the innermost terms in the Fourier Transform 
Decomposition: 

n«0 

10 Each such term has the form of a Discrete Fourier 
Transform. Therefore the efficient Fast Fourier Transform 
(FFI) algorithms can be applied to each of these terms. Note 
that there is an additional term in the exponent, (Tpk/Tpi), 
which is integral to the invention. But this is a constant for 

15 the series and does not change the FFT structure. 

Another point to consider is that although the outer 
summation appears to be unbounded, it is bounded in any 
physical situation. That is, only a finite number of prime 
scales need to be included in the D Scale, because no 

20 physical process can have an infinite bandwidth. Practical 
constraints also dictate realistic bandwidths as well. Thus, a 
voice application does not need generally need to consider 
a base bandwidth of higher than 20 KHz. 
Yet another point is that this equation quantifies just how 

25 much is thrown out by the standard approximation of a 
Fourier Transform by a Discrete Fourier transform. Indeed, 
all but one prime scale is disregarded in the usual approxi- 
mations. Therefore, the noticeable accuracy that must result 
from including the multi-resolutions will have real impact 

30 on the many commercial applications that are based on the 
Fourier Transform. 

Finally, because this decomposition is derived from first 
principles, it preserves the phase as well as amplitude of the 
original Fourier transform. 

35 There are many practical ways to sample data using the D 
Scale, and then to calculate the Fourier transform. The 
varied approaches depend on operating circumstances, such 
as real-time vs. post-data acquisition analysis, or available 
a-priori knowledge or not, as well as other cases. An 

40 embodiment is given below to highlight one manner by 
which the D Scale is used to calculate the Fourier Transform 
of sampled data during signal processing. 
Practical Example of Processing a One Dimensional Signal 
An example is presented that implements the use of the D 

45 scale in signal processing. Although the invention is appli- 
cable to multidimensional signals and images, etc. from 
many fields, a preferred embodiment will be described for a 
weak, sudden, signal over a telephone wire pair having a 
small bandwidth. 

50 The typical signal of this type has an exponential decay. 
This particular shape has a spectrum that is directly deriv- 
able from its Fourier transform. Its exact spectrum serves as 
a benchmark to compare the relative accuracy of the usual, 
textbook, Discrete Fourier Transform approximation against 

55 this invention. 

Two results emerge from this example. One is that the 
usual, textbook, approach requires many more uniformly 
spaced samples, yet still does not generate an accurate 
replica of the benchmark spectrum of the sound signal. 

60 Second, this invention uses far fewer points and does yield 
an accurate spectrum, comparable to the original. 

This benchmark then validates this invention's method- 
ology and decomposition equation. 

FIG. 9 shows the signal whose spectrum will be extracted. 

65 In order to demonstrate the effectiveness of the invention the 
benchmark signal was reproduced in a closed analytic form 
using Mathematica. 
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FIGS. lOa-fc are the mathematically exact spectrum of 
FIG. 9 to which will be compared the methodology of the 
prior art and of the invention. 
Result of Prior Art Approximation of Signal 

The Fourier Transform cannot be generally applied 
directly to data sets, uniform or not. One barrier is that data 
sets cannot, in the vast number of cases, be expressed in a 
mathematically tractable form. More important, the expo- 
nential term in the Fourier Transform causes partial to 
complete cancellation of terms generated from uniform 
sampling, leading to inaccurate spectrum calculations. 

These limitations notwithstanding, the usual approach is 
to approximate the Fourier Transform by a Discrete Fourier 
Transform (DFT) on uniformly spaced data. However, if the 
uniform rate is lower than the maximum possible 
bandwidth, then vital information could be lost in the 
calculated spectrum. Therefore the uniform sampling rate is 
set to the Nyquist bandwidth. This employs only one term of 
many that the invention will utilize: 



F(btml7 pi )~-n £ h(nf p k )e- i2Knm l T ^ 



values. This step effectively mixes the values from the prime 
subscales. Then identical entries are removed. This occurs 
because each unit interval is a part of each scale, like 1,2, etc. 
Finally, as a Mathematica artifact, the last entry in the 
5 temporary, intermediate, scale is appended. 

Do[sp[i]»Nrtable[(n-l)/p[[i]], {n, 1, p[[i]]pt[[i]]}]], {i, 
Lcngth[p]}]; 

scaleTmp=Sort[Flatten[Join[Table[sp[i], {i, 1, length[p] 
io }]]]]; 

scale-Table[Null, {0}]; 

Do[If[UnsameQ[scaleTmp[[i+l ]], scaleTmp[[i]]]scale= 
Append[scale, scaleTmp[[i]]], Null], {i, 1, Length 
[scaleTmp]-!}]; 
15 scale-Append]scale, scaleTmp[[Length[scaleTmp]]]]; 

creates the D scale. The intervals or widths between adjacent 
D scale points is calculated: 



20 



The calculated spectrum is arrived at by using Mathemati- 25 
ca's Fourier function. As may be seen from FIG. 11, this 
differs markedly from the exact spectrum in FIG. 10. Note, 
in particular, that the amplitude at 0 is much greater than the 
expected value of 2. Also note that the curvature in the lower 
left hand corner differs from the exact spectrum. 
Description of the Non-Uniform Sampling Device 

A measuring device is used to generate real-time, non- 
uniformly spaced samples, of a sudden sound over a tele- 
phone wire pair. The basic wire-pair telephony audio range 
has a 4 KHz baseband bandwidth. But the sound is weak 
enough that the effective base bandwidth is about 6 Hz. 
There is no loss in generalization, using this low frequency, 

Non-uniformly spaced samples -are generated by a mea- 
suring device that dynamically changes its sampling rate 
based on the volatility of its most recent samples. That is, it 
generates more Samples when recent samples suggest higher 
volatility, and far fewer samples where the time series is 
relatively constant. It does this by superimposing initial 
uniform sampling rates with new sampling rates that it 
activates when volatility increases. This provides a multi- 
resolution sampling effect. FIG. 4 shows how this device can 
be constructed using off-the-shelf components. 
Invention's Spectral Determination of Non-Uniformly 
Sampled Input Signal 

Consider the original signal as shown in FIG. 9. Its most 
volatile region is between 0-2 followed by 2-5 time units. 
Lastly, the region between 5-11 is relatively constant, with- 
out much information, therefore the non-uniform sampling 
device constructed above is programmed to partition the 
incoming data into three D scale subscales, at rates 1/11, 1/7, 
1/3. 

The initial rate is 1/3 for the entire 11 time units. The 1/7 
rate is applied only until the volatility reaches a low enough 
threshold, or 5 time units. Finally, the highest sampling rate 
of 1/11 is only applied for the first 2 time units, where 60 
volatility is highest. 

The D scale is constructed from the selected prime 
sampling rates, noted in the previous section. The first step 
creates the D scale's prime subscale vectors, in this case 
three. The Mathematica operation calculates decimal 65 
values from the fractions. These subscales are then joined 
and then most importantly, sorted from low to high numeric 



dclta-Tablct 5 calc[[i+l]]-scalcl[ill {*, 1, Length[scalc}-1}]; 

The values of the original signal are calculated at the D 
scale points. Non uniform sampling at the scale points is 
determined by 

/t=-Bible[Exp (-.5 scale[[i]]] UnitSte P [[*]]J {/, 1, Length! scale]- 



The D scale is indexed into its prime subscales. This is a 
programming artifact, not intrinsic to implementing the 
30 algorithm. The aim here is to effectively unmix the values in 
the D scale (due to the sort operation above), and place their 
D scale locations, or indexes, into their associated prime 
subscale. Each subscale vector is generated by first locating 
its members in the D scale. A member is identified if the 
35 product of this subscale's prime value and a scale value is an 
integer. Then, that member's D scale index is placed in the 
subscale array by: 



Clear[b] 
Do[ 
Do[ 

14 

IntegerQ[ 
Rationalize[ 

scale[[n]] p[[m]] 

] 
1 

b[m]=Append[b[m], n] 
Null 

], {n, Length[ scale]} 
], {m, Length[p]} 



45 



50 



55 



]; 

This step adjusts the index vectors created above. It 
removes those boundary scale entries like l,2,etc. that were 
placed in each subscale. These scale points are assigned 
arbitrarily to one of the subscales only, in this case the first. 
Then it is removed from the other subscales. Again, this is 
a programming artifact and would not be necessary if 
another implementation approach were taken. 

Finally, the implementation of £ S is: 
Do[f[l]« 

Table[ 

{2Pim/pt[[l]L 
Abs[ 
Sum[ 
Sum[ 
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h[[ c[k][[n]] ]] DFT approximation is really one of many DFT components 

delta[[ c[k][[n]] ]] " of the Fourier Transform above. That is, the invention shows 

Exp[ that the standard approximation throws out most of the 

-2Pi (pt[[k]]/pt[[l]]) relevant terms. 

(nm)/(p[[k]]pt[[k]]) I s g- The invention is faster because fewer points are used in 

] the calculation. And the invention shows how the calculation 

, {n, 1, Length[c]]-1} is comprised of decoupled terms, each with the FFT form. 

] b Therefore the FFT efficiency can be applied, as well as 

] parallel processing optimizations. 

} h. The invention also enables the phase of the spectrum to 



] 



]; 



f m 0 LenfithTcrin) ^ e accuraleiv calculated, because it relies on first principles. 

This is quite important for many applications, where phase 
(I 1 Len thf 1) is the determining factor in designed or observed behavior. 

LPl/ - Although the invention has been described in view of 

particular embodiments, the scope of the invention should 



the decomposition equation is implemented, now that the is be determille d by the fo'llowing claims: 

data structures are in place. Following the structure of the What is claimed is* 

implementation, a set of frequency amplitude vectors is L A method for procmdng a samp ied copy of a physical 

generated, from the sets of sampled values in the D scale's signal for signal proccssingj primarily digital signal 

prime subscales, again in this case three. The frequency processing by 

vector is denoted til], where 1 is indexed from 1 to Length 20 lin g the si , „ discrele mls of a D sca , 

[p], or time in this case. The D scale values are fetched wherein sM „ ^ ^ ^ um0Q of ^ of 

S ,t C[] , Ve 7^ e ^.^fP^^^P Th^' points, each set being points separated by regular 

b[ c[kl[n]] ] is the value of the original signal 'h at the scale m , ervals of y where is a ^ 

pomt located by its location m prune subscale vector c[k], at 2 A method for ali lhe Fourier transform of a 

position n in that vector. * . • 1 ■ i • • 

v _ , f . 1 . 25 physical signal comprising 

The amplitudes of these newly generated frequency vec- : . ,. t , , . . . . A ,. 4 . _ 

tors 51] are plotted, to highlight the underlying approach of « ^P' 11 * the P h ^ lc n al "t-duwete pomte of a D 

generating each subscale of the spectrum's D ULle. The ^ f where ' n siud . D com P nses "* umon ° k n 

Mathematica 'N' operation calculates decimal values from 56,5 ° f ?° ints \ ea f cl ? / set b ° mg P olnts se P arated u b y 

the complex numbers whose real and imaginary components 30 *S ular lnlervaU of 1/p « where p « 15 a prune number ' 

are expressed as fractions. 

Finally, the spectrum is plotted by juxtaposing the con- 00 transforming said sampled values into points in 

tributions from each of the spectrum D scale's prime sub- another parameter domain by a series of Discrete 

scales (each of which is shown in FIGS. 12-14). , ™ uner | ra ° sf °™?- 

Note that it accurately replicates the original spectrum 35 „ 3 ^ meth ° d °5 C \*™ 2 > wherein Said DlScrete Fourier 

from 0-10. Specifically, it has a value of 2 at 0, and its J'^™* m calculated usm % Fast Founer Tr ansforms 

curvature is very similar to the original. Thus, the value of ' V^t , * r , ^ 

0.5 is reached at about 1.75, as in the original. And the value u 4 ' method for generating the Founer transform of a 

of 1 is reached at about 1 on the x axis, again as in the P hvsical sl S nal of claim 2 > further comprising 

original spectrum. See FIG, 15. Thus the result of the present 40 ( a ) ^P^g me difference between the physical signal 

invention is excellent compared to the prior art. and the sampled, transformed result, 

Finally, the spectrum was recalculated (see FIG. 16) using ( b ) determining from the difference computed in (a) 

more data samples, at higher rates. The primes used were: whether to increase the resolution by recursing and 

p={29, 19, 11} instead of p={ll, 7, 3} increasing the set of distinct primes. 

The sampling time intervals remained the same as in the 45 5 - ^ method of claim 2 wherein kernel functions are 

prior case. Note that it also generates an authentic replica of integrated in the transform. 

the original spectrum. But, as expected, it is smoother than 6 - ^ method of claim 5, wherein the kernel functions are 

the previously generated one because more sample points chosen to produce the Radon Transform, Hartley transform, 

are applied. The point at value 0 is not indicated, although ^d Discrete Cosine Transform. 

it would be the correct value of 2, because the code used so 7 - ^ method of claim 6 wherein said method is imple- 

inadvertently started from the next to first index, not the first mented m a re al-ume data measurement environment, 

scale index. In practice the difference between a physical 8. A method for signal processing, primarily digital signal 

signal and the sampled transformed result may be calculated processing, by sampling the signal at discrete points and 

and from the difference it may be determined whether to transforming the values into an N dimensional Fourier 

recurse the sampling process while increasing the set of 55 transform, comprising 

distinct primes. selecting discrete points of a D scale at which sampling 

Advantages of the D Scale occurs, wherein said D scale comprises the union of n 

The advantages of using the Druck Scale are apparent se te of points, each set being points separated by 

when comparing the benchmark, textbook and invention regular intervals of l/p„ where p„ is a prime number 

calculation results. 60 transforming the sampled values into the coefficients of a 

e. The number of sample points required is dramatically Fourier transform of the signal, 

smaller than using the textbook approach. For it required a determining the accuracy of the transform and if neces- 

high data sampling rate, at twice the maximum bandwidth, sary recursing and increases the set of distinct primes, 

throughout the measurement period. 9. The method of claim 8 wherein said transformation 

f. The results are much more accurate because it does not 65 employs a Fast Fourier Transform analysis of the data, 
rely on the coarse Discrete Fourier Transform approxima- 10. The method of claim 8 wherein said method is 
tion. In fact, as the closed form expression indicates, the implemented by running a series of clocks ticking at inter- 
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vals characterized by l/p M for the n'* clock, and sampling the 23. The method of claim 8, wherein the method is applied 

signal at each clock tick. to data acquisition in one or more dimensions. 

11. The method of claim 10 wherein the signal is char- 2 4. The method of claim 23, wherein the method is 
actenzed by multiple parameters, the steps of the process .. , , . , , . , 

, . ; j t i j*j j *i r i_ applied after data has been acquired. 

being repeated separately and independently for each param- c *i _ . . 

eter 25. The method of claim 24 wherein the acquired data is 

12. The method of claim 11 wherein a sampling space is produced by NMR or tomography devices. 

formed as the Cartesian product of one dimensional spaces. 26. The method of claim 24, wherein the acquired data is 

13. The method of claim 10 wherein the signal is char- produced by seismographic measuring devices, 
acterized by multiple parameters, the steps of the process 27. The method of claim 8, wherein the method is applied 
being repeated separately and concurrently for each param- i n real-time. 

eter * 28. The method of claim 8 wherein the D scale is used to 

14. The method of claim 13 wherein a sampling space is ^on or creat6 mm h of physical> ; u . 
formed as the Cartesian product of one dimensional spaces. ^ j atterned data 

15. The method of claim 8 wherein the physical signal is -« ~- \ / * . 

characterized by multiple parameters, the steps of the pro- « 29 - ™ c method of claim 28 > wherein said P artltl0n 

cess being repeated separately and independently for each contains periodically patterned data from physical samples, 

parameter. 30. The method of claim 29, wherein said data is band- 

16. The method of claim 15 wherein a sampling space is limited data from the physical samples. 

formed as the Cartesian product of one dimensional spaces. 31. The method of claim 30, wherein each partition has a 

17. The method of claim 8 wherein the signal is charac- 20 different resolution such that the partitions form a nested set 
terized by multiple parameters, the steps of the process being Q f resolutions 

repeated separately and concurrently for each parameter. 32 ^ method of ckim 28, wherein each partition has a 

18 Tlie method of claim 17 wherein a sampling space is resolution such that the rtitions form a nested 

formed as the Cartesian product of one dimensional spaces. . ^ c i ■ . 

19. The method of claim 8 wherein said signal is a time 25 set of non-overlapping points. 

series originating an electrical, optical, or other physical . 33 ' ^ mcthod of claim 32, wherein each partition has a 

process, wherein the sampling occurs at time or other different resolution such that the partitions form an infinite 

parameter values selected from a D scale. " nested covering set of non-overlapping points. 

20. The method of claim 8 wherein said signal is a series 34. The method in claim 28 wherein a maximum and 
resulting from data processing on any kind of physical data, 3Q minimum resolution of the nested partitions is pre- 
wherein the sampling occurs at time or other parameter determined such that neighboring points in the combined 
va h es Jf lecled fr ? m , a , D partitions have an interval between the specified minimum 

21 The method of claim 8 wherein the sampled object is and maximum resolutions. 

a multidimensional object, and wherein the signal comprises ~, mo ,u rt/1 rt f e „ r . ■ Antn „ 

N independent variables each of which is measured at D ™ e meth °? ° f claun 8 ; where * the u data f*P™?* 

scale values comprising a multi-dimensional coordinate 35 random noise, and the n sets of points filter the random noise 

system. mt0 periodic components. 

22. The method of claim 21 wherein the multidimensional 36 The method of claim 35 further comprising the step of 

object is embedded in a multidimensional space and the D processing the periodic components, 
scale values parametrize the space in which- the object is 

embedded. ***** 
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